Hydrodynamic orienting of asymmetric microobjects under gravity 
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It is shown that nonsymmetric microobjects orient while settling under gravity in a viscous fluid. 
To analyze this process, a simple shape is chosen: a non-deformable 'chain'. The chain consists of 
two straight arms, made of touching solid spheres. In the absence of external torques, the spheres are 
free to spin along the arms. The motion of the chain is evaluated by solving the Stokes equations with 
the use of the multipole method. It is demonstrated that the spinning beads speed up sedimentation 
by a small amount, and increase the orientation rate significantly in comparison to the corresponding 
rigid chain. It is shown that chains orient towards the V-shaped stable stationary configuration. In 
contrast, rods and star-shaped microobjects do not rotate. The hydrodynamic orienting is relevant 
for efficient swimming of non-symmetric microobjects, and for sedimenting suspensions. 

PACS numbers: 47.20.-k, 82.70.-y 



I. INTRODUCTION 

In many biological, medical and industrial applica- 
tions, it is of interest to predict theoretically what is the 
sedimentation velocity of small conglomerates of micro- 
particles under gravity in a fluid 0,1, a ii, and how 

the settling speed can be enhanced or decreased, by a 
suitable modification of the configuration, or directly by 
a change of the relative motion of the particles. This issue 
is especially important for mechanisms of effective swim- 
ming, recently intensively investigated for biological sys- 
tems as well as for artificial micro- and nano-swimmers. 

Swimming patterns of various microorganisms have 
been extensively investigated experimentally and theo- 
retically [i, 0, i, i, M, M, m, E [ll- Microorganisms 
propel themselves owing to a periodic change of their 
shape and possibly also its orientation in space. Often the 
core cell does not deform, and the shape is changed owing 
to waving, ondulating or rotating flagella. Typical sizes 
of bacteria, spermatozoa or algae lie in the range from 
1 to 200 /im, and their swimming speeds are usually up 
to several hundred /im/s. For such microobjects moving 
in aqueous environments, the fluid inertia and the Brow- 
nian motion are irrelevant (111] . Therefore, a theoretical 
model of swimming should be based on hydrodynamic 
interactions 15] between individual parts of the microob- 
ject, following from the stationary Stokes equations and 
the appropriate boundary conditions at the surface of the 
swimmer. Typically, in the swimming problem a periodic 
sequence of flagella shapes and the corresponding trans- 
lational and rotational velocities of their parts relative 
to the core cell are known as functions of time, as well 
as the total force and torque exerted on the cell center. 
The task is to determine the translational velocity of the 
center, and also the cell's angular velocity. There exist a 
number of models of freely moving swimmers (the total 



force and torque vanish). 

However, the microorganisms are often denser than the 
water in which they swim, by a few percent for the algae, 
approximately 10% for bacteria and 30% for spermato- 
zoa, and the mass distribution can be non-uniform pd| . 
The gravitational force is essential for explanation of ori- 
entational mechanisms (gravitaxis) observed experimen- 
tally e.g. for algae [l6l.ll7j. In general, hydrodynamic in- 
teractions would tend to orient non-symmetric microob- 
jects settling under gravity. This effect, certainly impor- 
tant for swimming, will be investigated in this paper. 

We focus on a very simple model: a non-symmetric 
'chain' of three identical spheres, with two pairs at con- 
tact, but not the third one. Owing to the lubrica- 
tion forces, the shape of the conglomerate is fixed. In 
Sec. ini the accurate spherical multipole method [H, [l^ 
of solving the Stokes equations is introduced, with the 
lubrication correction for the relative motion of close 
surfaces [20|, and the HYDROMULTIPOLE numerical 
code [13] ■ In Sec. IIIIl we use this method to evaluate 
the translational, rotational and spinning velocities of the 
microobject, determine how it orients while settling un- 
der gravity, and find the stable stationary configuration. 
In Sec. IIV| we study how the settling speed depends on 
shape, by comparing sedimentation velocities of all the 
stationary configurations of three spheres [2l|, [S^, H^. 
We also investigate how accurate is the point-particle 
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approximation 2M- In Sec. |V] we summarize the results 
obtained for the chain made of three spheres. We also 
check if chains with two arms made of a larger number of 
beads orient hydrodynamically, qualitatively in the same 
way as the simple three-sphere model. 



II. THE MODEL OF A MOVING 
ASYMMETRIC MICROOBJECT 

Consider three identical spheres falling under gravita- 
tional forces Fq in an infinite fluid of viscosity 77. The low 
Reynolds number is assumed for the corresponding fluid 
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flow. The fluid velocity v(r) and pres sure p{r) satisfy 
the stationary Stokes equations [isl. l25j|. 



TjV^v -Vp = 0, 



0, 



(1) 



with the stick boundary conditions at the surfaces of the 
spheres and no fluid flow at inflnity. Therefore the trans- 
lational Ui and rotational fli velocities of each sphere 
i = 1,2,3 are linear functions of the force Fq, 
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The 3x3 mobility matrices /i*^ and are to be found 
as functions of the relative positions f / — rj of the sphere 
centers. 

In the following, as in Ref. ^25|, particle positions f ; 
will be normalized by the sphere diameter d, translational 
velocities Ui by the Stokes velocity. 



Us = |Fo|/(3^?7rf), 



(4) 



angular velocities fli by 2Us/d, and time t by two Stokes 
times, 2rs, with 



TS 



d/{2Us). 

The corresponding dimensionless quantities are 
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and the dimensionless Eqs. Q-® read, 
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(13) 
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In our model it is assumed that each sphere touches 
another one. Once it happens, the spheres remain at 
contact owing to lubrication forces [26[. Such conflgu- 
rations will be called 'chains'. As illustrated in Fig. [U 
positions of the sphere centers are parametrized by the 
angle a between the chain links. 



(0,0,0), 
{x/2,0,z), 



(15) 
(16) 
(17) 




FIG. 1: Parametrization of a chain configuration. 



with X ^ 2 sin(Q:/2) and z ~ cos(q!/2). Orientation of the 
gravitational force is arbitrary, Fq = {Fq^, Foy^ Fqz)- 

Evaluation of the mobility matrices of chains requires 
a special treatment, because now there are at least two 
contact points. Moreover, it is essential to specify more 
precisely what is meant by "the contact" . The flrst pos- 
sibility is to assume that there is no external forces other 
than gravity exerted on each of the spheres. In this case, 
lubrication does not allow for any relative motion of the 
touching spheres except their spinning along the line of 
centers ^] . In the following, such a contact will be called 
"beads". The second option is to "glue" the touching 
spheres, imposing on them external torques, which pre- 
vent them from any relative motion. Such a contact will 
be called "rigid" . In Appendix [Aj it is explained in de- 
tails what is the difference between the mobility matrices 
for both types of the chains. 

The essential task of this paper is to determine the 3x3 
mobility matrices /x*^ and fi^j: both for chains made of 
beads and the rigid ones. Because of hydrodynamic in- 
teractions of the close surfaces, we have to go beyond the 
point-particle approximation, and even beyond superpo- 
sition of two-body mobilities [l^ [l^ . 

Therefore we evaluate the three-sphere mobilities nu- 
merically by the multipole expansion [l^ [l^ . The algo- 
rithm from Ref. [l^ and its accurate numerical FOR- 
TRAN implementation HYDROMULTIPOLE ^ are 
applied. The accuracy is controlled by a varied order of 
the multipole truncation L (see Ref. [HI for the deflni- 
tion of L and Refs. [H, H^l for discussion of the accuracy 
estimates). In this paper, we take a large value L = 6. 
The results will be presented in the next section. 



III. MOTION OF A CHAIN 

It is convenient to describe the motion of a chain (rigid 
or made of beads) referring to its center of mass, R — 
(ri + r2 + T3)/3, because in this case the total external 
torque vanishes. The task is to flnd its translational and 
rotational velocities, 

U = {Ui + U2 + U3)/3, (18) 

n = n2^{ni + n2 + n3)/3. (i9) 
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For the chain made of beads, the spinning speed u/"'"-''-^ 
also needs to be determined; from the symmetry of the 
system it follows that 



— Vti — Q,2, and 



heads 



(20) 



with the unit vector rij = rij 



j I , and the standard 
notation for the relative positions, r^j = — r^- . 

The translational and angular velocities of the chain 
depend linearly on the total gravitational force acting on 
the chain, F = 3Fq, 
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with the mobility matrices of the chain to be found. Ow- 
ing to the symmetry, in the frame of reference shown in 

Fig.m 



Aii(a) 
/i2(a) 
^l3ia)^ 

Mb(a) 
-pia(a) 











-( M-(a) ) 



(24) 

(25) 
(26) 



In the above equations, the factor 1/3 has been intro- 
duced. With this choice, in the frame of reference shown 
in Fig. [U and with the adopted normalization 
the mobility coefficients for the chain are just equal to the 
corresponding velocity components, ^ii = Ux, = C^i 



= Uz, fJ^a = —^y and jib — ^x- The physical meaning 
of the coefficients is indicated in Fig. [5J 

The mobility coefficients in Eqs. p4)) - ((25)) are deter- 
mined from combinations of the three-sphere friction co- 
efficients for the individual spheres, as outlined in Ap- 
pendixIXl In general, the mobility coefficients for a chain 
of beads differ from those for a rigid chain at the same 
configuration and in this case they will be denoted by 
the corresponding superscripts. Thus, fJ^'^'^'^^ ^ /i"^*'' 
and ^ /i"^*''. For a rigid chain, there is no spin- 

ning and /i"^*'' — Lu"s^'^ = 0, while the chain of beads 
does spin, ^^J"-'^'^ — ^beacU _i_ Q However, for a rigid chain 
and the chain of beads at the same configuration, /ii, 
/i3 and [la are identical and those coefficients will not be 
marked by any superscripts. 

The translational motion will now be determined. The 
corresponding mobility coefficients are evaluated numeri- 
cally and plotted in Fig.[3]as functions of the angle a. For 
a ~ /i, the sphere centers form the equilateral triangle 
with three contact points; in the following, this configura- 
tion will be called 'a star'. For a = tt, the sphere centers 
are aligned; this configuration will be called 'a rod'. For 
112 and /Z3, when the base of the chain is perpendicular 




FIG. 2: Degrees of freedom of chains at the characteristic 
orientations with respect to gravity. Top: gravity along x; 
middle: gravity along y; down: gravity along z. 
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FIG. 3: The chain translational velocity components. 



to gravity, then its settling velocity is a decreasing func- 
tion of the apex angle a, in agreement with the intuitive 
prediction that stretching the arms should increase the 
friction force. For /ii, when the base of the chain is paral- 
lel to gravity, it is intuitive to expect the opposite effect: 
with the increase of a, the chain aligns along gravity, 
and its resistance is weaker. This is indeed observed for 
a wide range of the larger angles, except those relatively 
close to 7r/3. 

For a given shape (fixed a) , we now compare the mag- 
nitude of /ii, /i2 and 113. Settling is always the fastest if 
gravity is along the base of the chain, and the slowest if 
gravity is perpendicular to the plane of the chain. Notice 
that in the last case, the chain made of spinning beads 
sediments faster than the corresponding rigid one. 
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The spinning speed is plotted in Fig. [4] It reaches a 
maximum for a very smaU value of a, and then slowly 
decreases. For a wide range of the angles between the 
chain arms, the spinning surfaces move with velocities 
which are still around 5% of the settling speed. 



0.12 




a 



FIG. 4; The spinning speed of the chain of beads. 

Next, we evaluate the chain rotation. The correspond- 
ing mobility coefficients are plotted in Fig. [5] as functions 
of a. Each component of the angular velocity is zero 
for stars and rods, and has a maximum at intermedi- 
ate values of the angle between the arms, smaller than 
2tt/3 for a rigid chain. Therefore the speed of the hy- 
drodynamic orienting is very sensitive to shape. It is re- 
markable that spinning of the beads has a profound effect 
on the chain rotation. The spinning increases the angu- 
lar velocity by at least a factor of two in comparison to 
the corresponding rigid chain at the same configuration. 
Moreover, we observe a much wider range of the chain 



0.08 




a 



FIG. 5: The chain angular velocity components. 



shapes for which the rotation is significant. Indeed, the 
maximum is shifted to larger values of a, above 27r/3, and 
there is a qualitative difference at small angles, where the 
lubrication interactions between the spinning beads keep 
the chain rotating. 

When the close surfaces of spheres 1 and 3 move 
with respect to each other, then the mobility coefficients 
^beads ^ ^beads ^^^^ ^beads (jgcreasc rapidly with the de- 
creasing a only if a — 7r/3 becomes extremely small, 
as seen in Figs. [31 [SI and [D This is the typical lubri- 
cation interaction of very close surfaces in relative mo- 
tion [2^. Actually, the relative mobility coefficients de- 
crease to zero as the inverse logarithm of the gap size. 
This scaling can be seen in Fig. [S] where ^^"■'^'^ , 
and ^'^"■d-^ are plotted as functions of [— l/ln(a — 7r/3)]. 
In fact, to account for the relative motion, we plot 
A*2(o) ~ fJ'2{'^/S) rather than /i2(o), with the zero-gap 
mobility /X2(7r/3) — 1.739 (see next section for the deriva- 
tion). Numerical results are not available when the size 

0.03 1 . 1 




-l/ln(a-7r/3) 



FIG. 6: The mobility coefficients for a spinning chain of beads 
with a small a. 

of the gap between two surfaces is comparable with the 
numerical accuracy. For smaller values of a, asymptotic 
expressions analogical to Eqs. (12)-(13) in Ref. [22] could 
be derived using the same procedure. Notice that in the 
lubrication regime, the spinning speed of the beads is sig- 
nificantly larger than the angular velocity of the chain. 

It remains to discuss the main issue of this paper, that 
is how the chains orient with time. Taking the frame 
of reference (x', y' , z'), in which Fq = (0,0,-1), and 
using the spherical coordinates, we parametrize the chain 
symmetry axis z (see Fig. [T] for the sense of z) by the 
angles 9 and (p. We now have cos 9 = —Fq ■ z, and we are 
interested in the time evolution of the angle 9. In general, 
it depends on the angle ip between the line of centers of 
spheres 1 and 3 and the unit vector eg = dz/d9. From 
Eqs. ([HI) and ^ we obtain, 

9 — —{fJ-a cos^ ip + ^J-b sin^ V) sin 9. (27) 
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For 7r/3 < a < TT, the calculated mobility coefficients fia 
and are positive, and therefore from Eq. (j27p it imme- 
diately follows that 9 evolves towards zero, e.g. towards 
the chain-axis antiparallel to gravity. This orientation of 
the chain will be called 'V-shaped'. It is the only stable 
orientation of the chain with 7r/3 < a < tt. 

Two examples of a one-dimensional dynamics are re- 
covered from Eq. (j27p for two symmetric cases with 
■0 = and V' = 7r/2, respectively. For V' = 0, Eq. ([TT]) 
reads 9 = — /ia sin 0. This dynamics is illustrated in 
Fig. [21 with its upper C-shaped configuration, which cor- 
responds to 9{t = 0) = tt/2, evolving towards the bot- 
tom one (V-shaped). The sense of the rotation is il- 
lustrated by the arrow. For if; — tt/2, Eq. (P?]) reads 
9 ~ —jib sin 9. This dynamics is also illustrated in Fig. [21 
now with its middle plane configuration, which corre- 
sponds to 9{t = 0) = it/ 2, evolving towards the bottom 
one (V-shaped). The sense of the rotation is illustrated 
by arrows of different lengths. 

To complete the analysis, we still need to compare a 
typical time scale of the hydrodynamic orienting to that 
characteristic for the settling motion. We therefore find 
the angles aa,max and ab^max which correspond to the 
maxima of fj,a and fj,b, respectively, and then calculate 
the ratio of both time scales, 27r/ii//ia and 27r/i2//^6, at 
aa,max and ab,inax, respectively. For the rigid chain, both 
ratios are of the order of 300. For the spinning beads, 
27r/i2/A'6 and 27r/i2//iuj at Umax are of the order of 150. 
Therefore the characteristic time scale of the hydrody- 
namic orienting is at least two orders of magnitude larger 
than that of the gravitational settling. Reorientation of 
the sedimenting nonsymmetric particles is significant for 
such systems which stay under gravity for a sufhciently 
long time. 



IV. STATIONARY CONFIGURATIONS 

The goal of this section is to show what are the chain 
configurations, which do not orient under gravity. Set- 
tling speeds of such stationary configurations will be in 
addition compared to the translation velocities of other 
stationary configurations of three spheres, with the em- 
phasis on those with the spinning particles. 

By definition, at a stationary configuration the spheres 
have equal translational velocities, 

t/^ = U. (28) 

Such a configuration is an equilibrium solution of the 
dynamics of the relative positions. Notice that U is time- 
independent. In our case, obviously, i = 1,2,3. 



A. Stationary chains 

For chains, the equilibrium condition ([28)1 is equivalent 
to the relation, 

rj 0, (29) 

which takes the form 

Ha{a) = ^lb{a) Fy ^ 0, (30) 

if Eqs. and ([^5)1 are applied. According to the nu- 
merical results plotted in Fig. [51 Eq. (|5D)) has the solu- 
tions, 

(i) an arbitrary a, and F^ — Fy — 0, 

(ii) an arbitrary F, and a = 7r/3, or tt. 

The condition (i) corresponds to the vertical chain equi- 
libria, found in Ref. [22| and sketched below in Fig. [71 
These are the chains with the symmetry axis parallel to 




FIG. 7: Vertical chain configurations. 

gravity. The condition (ii) corresponds to stars and rods, 
sketched in Fig. [H 




FIG. 8: Star-shaped configurations with a = n/3 (left) and 
rods with a — n (right). The direction of gravity is arbitrary. 

Notice that from the symmetry it follows that the 
beads of the stationary chains do not spin, w = 0. The 
settling velocities U of the stationary chains will now 
be evaluated. The results will be also compared with 
the point-particle model. For touching spheres, such an 
approximation has to take into account additional con- 
straint forces, which do not allow the points, which ap- 
proximate the touching spheres, to change the interpar- 
ticle distance [2^1 • 

Consider first the vertical chains. From symmetry with 
respect to reflections x,X —>■ ~x,—X and y — > —y it 
follows that U is vertical. Values of the settling velocities, 
U{a) = ^3 (a), evaluated in Ref. [l^l, and replotted here 
in Fig. [21 span the range /i3(7r) <U< /i3(7r/3), with 

^3(7r) = 1.63045819, (31) 
/X3(7r/3) = 1.9022670. (32) 

Similar values follow from the point-particle approxima- 
tion with constraints, 29/18 < < 229/112. 
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rod ± F 

ceo 


vertical chain 


star II F 

QD 


star ± F 

CD 


kissing 

oo o , 
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ring 

oo 


spheres 


1.95 


1.63 


1.63-1.90 


1.90 


1.74 


1.75 


1.74-1.79-1 


point-particle 
approximation 


2.00 


1.61 


1.61-2.04 


2.04 


1.75 


1.81 


1.75-1.74-1 



TABLE L Settling velocities U of stationary configurations of three spheres. 



The settling velocities J7|[ and U± of the rods parallel 
and perpendicular to gravity are now evaluated for the 
subsequent values of the multipole order L < 30, and 
extrapolated to L — > oo. Then, 

U\\=fii{n) = 1.946299144, (33) 
Ui_ = M2(7r) = 1.63045819. (34) 

These values are again well approximated by the point- 
particle approximation with constraints, with = 

2 and ^p"™*" = 29/18. Both parallel and perpendicular 
rods settle down vertically, i.e. along gravity. The calcu- 
lated velocities agree with the previous experi- 
mental and numerical results P, Q . 

In a similar way we calculate velocities of the stars. 
Notice that owing to the symmetry with respect to rota- 
tion by 7r/3, the stars located in the vertical plane settle 
with the same velocity, independently of their orienta- 
tion. However, their settling velocity U\\ is larger than 
that of the stars oriented horizontally, U± . In both cases, 
the stars translate vertically. We evaluate, 

= ^i(7r/3) ^ 1.90226703, (35) 
[/^=^2(7r/3) = 1.73941260. (36) 

In the point-particle approximation with constraints, 
^po^nts ^ 229/112 and t/f = 7/4. The results ([23- 

(PS)) obtained for the stars improve the accuracy of the 
previous simulations [1], and agree well with the mea- 
surements in. 

In general, the stars and rods are inclined at a certain 
angle with respect to gravity. In this case, their veloc- 
ities are not vertical, and the components follow from 
Eq. ([M)) . A special case of such inclined stars was dis- 
cussed in Ref. [22] , where it was indicated by the dotted 
line in Fig. 15. In this "slanted equilateral chain" con- 
figuration, a line of centers was perpendicular to gravity. 

Using the dynamics derived in the previous section, We 
conclude that the only stable stationary solutions of the 
dynamics of chains are the V-shaped vertical chains. The 
hat-shaped vertical chains are unstable. Rods are neu- 
trally stable. So as the stars if the three contact points 
are kept. The stars are unstable against perturbations 
which separate out a pair of the touching spheres. 



B. Comparison with other stationary 
configurations 

The settling speeds of stationary chains will be now 
compared with the motion of other equilibrium configu- 
rations of three spheres. At the equilibrium, the triangle 
formed by the sphere centers has the following shape, size 
and orientation with respect to gravity. 

• "Vertical chain" (a vertical isosceles triangle with 
the symmetry axis along gravity; the apex sphere 
touches each of the base spheres). 

• "Rod" (a straight line of an arbitrary orientation 
with respect to gravity; there are two contact points 
between the sphere surfaces). 

• "Star" (an equilateral triangle at an arbitrary ori- 
entation with respect to gravity; there are three 
contact points between the sphere surfaces). 

• "Kissing" (an isosceles triangle with the symmetry 
axis along gravity and the touching base particles; 
the distance between the contact point and the cen- 
ter of the apex particle equals 1.578634 diameter, 
see Ref. dij and [l^). 

• "Ring" (an equilateral triangle of an arbitrary side 
length, in the plane perpendicular to gravity, see 
Refs. [gllll). 

In Table HI the stationary configurations are sketched 
and values of their vertical velocities are listed, together 
with their approximation by point-particles with con- 
straints. We have demonstrated that for a small number 
of particles, the settling velocities of their stationary con- 
figurations can be within a few percent approximated by 
the point particle model with constraints. Notice that 
all the equilibria except the ring arc unstable, if arbi- 
trary perturbations are allowed, including separation of 
the touching surfaces [23 |. 

The ring is the only equilibrium configuration with the 
rotation of the individual spheres. It is therefore inter- 
esting to investigate if the spinning increases the settling 
velocity, as it has been observed for chains made of beads. 
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This problem will be discussed in details in a separate 
section. 



C. Stationary configurations with spinning 

We now focus on the stationary configurations called 
rings. The sphere centers form a horizontal equilateral 
triangle with an arbitrary length ^ > 1 of its side. Set- 
tling velocity U is of course vertical. Its value U = \U\ 
is evaluated and plotted in Fig. [H] as a function of £. In 
general, the particles are separated from each other; they 
touch only in the limiting case of ^ = 1, when the ring 
becomes the horizontal star, with U = 1.73941260. 

It is interesting to observe that the ring's settling ve- 
locity has a maximum for a very small gap between the 
sphere surfaces. The maximum is well-visible at the in- 
set of Fig. [S) Bracketing this maximum by the standard 
golden section search [291], we evaluate the correspond- 
ing values of £mu and U{£mu)- The subsequent multi- 
pole orders 1 < L < 28 are used, and the results are 
next extrapolated to L ^ oo, as in Ref. 12^. We obtain, 
ij^t^j = 1.01128 and [/(^Afc/) = 1.79394. For larger values 
of the settling velocity decreases to zero with I ^ oo. 
In Tab. HI the ring velocities at £ = 0, Imu and oo have 
been indicated. 

To check if the existence of the maximum is related 
to the spinning, we also consider a rigid system of three 
spheres at the same configurations, but with the spinning 
eliminated owing to constraint external torques. Veloci- 
ties of such configurations are also plotted in Fig. [51 They 
are systematically smaller than the velocities of the ring. 
Indeed, this example also indicates that spinning speeds 
up the rate of settling by a small amount. 

The point-particle approximation, U — \ + 3/(4i?), is 
also depicted in Fig. [SI Notice that for rings, no con- 
straint forces nor torques are applied, since the rings are 
also stationary solutions of the point-particle dynamics. 
It is clear that the point particle approximation is much 
closer to the rigid dynamics than to the spinning system. 

Finally, we evaluate the spinning velocities of the 
spheres. Here fJi =ilo r32, ^2 — ^0 ^"13, and ^l^ — ^lo r2i, 
with fio, plotted in Fig.fTUlas a function of the side length 
i. The maximum of is reached at the small distance 
between the sphere centers, (mq = 1.0923791, but not 
as small as £mu- At the maximum, the spinning velocity 
^o{£Mn) = 0.215537174, is as much as 12% of the set- 
tling speed, U{£Mn) — 1.7519797. At the maximum of ?7, 
the spinning is slightly smaller, with floi^Mu) = 0.17704. 
These two maxima are shifted with respect to each other, 
because U becomes larger not only by an increase of 
ilo, but also by a decrease of the distance between the 
spheres. 

Notice that for the ring configuration of the separated 
spheres, the maximal spinning velocity is two times larger 
than for the horizontal chain of the touching beads. As a 
consequence, the increase of the settling velocity is also 
twice as large. 




1 2 3 4 5 



FIG. 9: Settling velocity U of the ring (equilateral horizon- 
tal triangle) versus its side length £. Equilibrium (solid line); 
rigid system with constraints (dashed line) , point-particle ap- 
proximation (point line) and the horizontal star (o). Inset: U 
as a function of [—1/ ln{£ — 1)] for very close particles. 
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FIG. 10: Spinning velocity flo versus the size £ of the ring. 

V. CONCLUSIONS 

The goal of this paper has been to construct and study 
a simple model of a chain-like asymmetric microobject of 
a fixed shape, settling under gravity in a viscous fluid. 
The motion of such a system has been evaluated from 
the multipole expansion of the Stokes equations. The 
main results are the following. 

It has been found that asymmetric microobjects orient 
hydrodynamically while settling under gravity. This ef- 
fect is not observed for axially symmetric objects such as 
rods, neither for regular shapes such as equilateral trian- 
gles, here called stars. However, chain-like conglomerates 
made of two identical straight arms in general orient to- 
wards a vertical "head down" equilibrium configuration; 
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that is, towards a V-shape. This process is relatively slow 
in comparison to the cluster sedimentation, but definitely 
not negligible. We have checked that the hydrodynamic 
orienting is observed for chains made of a central particle 
with the attached two identical straight arms made of a 
certain number of spinning or rigid beads, for example 
ten spheres at each arm. The central particle may be 
the same size as the other ones, or larger, e.g. 10 times 
larger. 

It has been shown that freely rotating particles in 
chain-like conglomerates can spin even if their surfaces 
touch each other. The spinning can be even faster than 
the chain rotation. The spinning particles speed up the 
conglomerate settling by a few percent, and they signifi- 
cantly enhance its tendency to orient vertically, in com- 
parison to the rigid body. Spinning also speeds up the 
settling of stationary configurations of spheres separated 
from each other. 

Naturally, the hydrodynamic orienting found in this 
work is important for efficient swimming of microorgan- 
isms which are more dense than the fluid. The results 
are also relevant for suspensions of chain-like conglom- 
erates settling under gravity. On the long time scale, 
while reorienting takes place, the suspension structure 
and settling speed may change, leading to ordering of the 
sediment and possible applications to segregation and fil- 
tration techniques. 



APPENDIX A: HOW TO EVALUATE MOBILITY 
OF A CONGLOMERATE? 

For the system of N particles separated from each 
other, the 6x6 mobility matrices /i.^ form a 6iV x 6A^ 
tensor, which is evaluated as the inverse of the 6A^ x 6A^ 
friction tensor, made of 6 x 6 friction matrices ^j^ . The 
latter relate the external forces and torques, Fi and Ti, 
exerted on a particle i = 1, N, to the translational and 
angular velocities, Uj and flj, of a particle j = 1, ...,N, 




(Al) 



In a conglomerate, the particles touch each other, and 
some of the ^.^ components become infinite. Instead of 
using Eq. (jAip . it is therefore necessary to eliminate the 
forbidden degrees of freedom (relative motions). It is 
done by constructing the conglomerate friction as the 
sum of the relevant combinations of Cij only. 



N N 



(A2) 



i=l j=l 



The 6x6 friction matrix ^ is finite at the contact, be- 
cause the combinations of Cij in Eq- (|A2p correspond to 
motions, which are free from the lubrication singularities. 
By inverting C, we obtain the conglomerate mobility. 



In the following, we construct the operators Pi for two 
different types of the conglomerates. First, a rigid system 
is considered, for which all relative motions are excluded. 
Then, a conglomerate with spinning particles is analyzed. 



1. Hydrodynamics of a rigid system 

A rigid motion of a conglomerate made of N spheres 
is characterized by the translational velocity U of an ar- 
bitrary chosen center of reference R and the conglom- 
erate rotational velocity ft. Then, the translational 
and rotational velocities of the individual sphere centers, 
i — 1, . . . , TV, are given as 

U, = CZ + n X (r, -R), (A4) 
= n, (A5) 

or equivalently. 





(A6) 



where the 6x6 matrices are given by the relation 

/ I ptr \ 



P. = 



I 



(A7) 
(A8) 

The total force F and torque T with respect to the cen- 
ter R, exerted externally on the conglomerate, have the 
form, 



N 



i=l 
N 

Et 

i=l 



T, + (r, -R)xF,], 



or equivalently. 




N 




(A9) 



(AlO) 



(All) 



where stands for the matrix transposition. By insert- 
ing Eqs. and ((X6|) into Eq. (jAlip . we relate the 
total external force and torque on the conglomerate to 
its translational and rotational velocities, 




= c 




(A12) 



(A3) 



with the conglomerate friction ^ given by Eq. (|A2p . Writ- 
ting the conglomerate velocities explicitly it terms of the 
corresponding components of the mobility fi = C~^j 



(A13) 



and choosing the center of mass as the reference center, 
we obtain T = 0, and we recover the relations (PT|) - (P^ . 
with the rigid-chain mobility matrices /x** and /x^*. 
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2. Hydrodynamics of a chain made of beads 



Let us now consider a chain of three beads. The sphere 
2 touches the other spheres, but the spheres 1 and 3 
are separated from each other and therefore are able to 
spin along fi2 and i'32, respectively. The motion of the 
chain is characterized by the translational velocity U of 
an arbitrary center of reference R, the chain rotational 
velocity fl, and also by the two spinning velocities wi 
and UJ3. Then, the translational and rotational velocities 
of individual sphere centers are given by the relation. 



U^ = U + ftx{r, 

O2 = o. 



1,2,3, 
1,3, 



(A14) 
(A15) 
(A16) 



or equivalently. 



= P^ 



1,2,3, (A17) 



with the 6x8 matrices Pi defined with the use of the same 
Eq. ((X8|) for Pf, but now differently than in Eq. (M\i . 



I P|'~ 
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3 = 




(A19) 
(A20) 



Then, Eq. (|A2p is used to evaluate the 8x8 chain friction 
matrix ^. It relates the the chain velocities U, fl, loi and 
UJ2, to the total external forces and torques (jA9P - (|A10p . 
and the torque components ti—Ti- ri2 and ts = • 7*32, 



(A21) 



Evaluating the chain mobility fi — C , we obtain, 



(F] 






T 






ti 




Ll>i 






\UJ3 J 







(F\ 


n 


= M • 


T 




h 


\UJ3 ) 







(A22) 



Pi = 



I Pf 
I ri2 



(A18) 



In our system, the spheres are identical, R is the center- 
of-mass position, T = and ti = — 0. Therefore, 
LOi = — wa = u) and we obtain the relations ([2T|) - ([23)) . 
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